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This paper discusses the scattering losses of dielectric slab waveguides 
that are caused by higher-order grating lobes of the sinusoidally distorted 
core-cladding interface. The results of this paper are used in a companion 
paper to evaluate the radiation losses of multimode guides with intentional 
mode coupling. An exact system of equations is derived for the amplitudes 
of all grating orders. This system is used to derive first- and second-order 
approximations that hold for small amplitudes of the sinusoidal interface 
distortion. The theory is used to derive formulas for the average power loss 
coefficient for first- and second-order scattering processes. 

I. INTRODUCTION 

Scattering losses in dielectric waveguide caused by core-cladding 
interface irregularities have been studied extensively by means of 
first-order perturbation theory. The principle result of this theory can 
be stated as follows: 1 Two modes with propagation constants jS„ and j8„ 
are coupled only if a Fourier component of the core-cladding interface 
function exists the mechanical frequency of which, <j>, satisfies the 
relation 

* = [ A ~ ft. I • (1) 

The propagation constants j3, and & may both belong to guided modes 
or one may belong to a guided mode while the other belongs to the 
continuum of radiation modes. Coupling of a guided mode to radiation 
modes results in power loss of the guided mode. This first-order coupling 
process is very strong and leads to high radiation losses if suitable 
Fourier components of the mechanical core-cladding interface irreg- 
ularity function exist. 

The result of first-order perturbation theory can be understood by 
viewing the core-cladding interface as a diffraction grating. 2 Since it 
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modulate? the phase of the incident wave passing through it, the dielec- 
tric interface acts as a phase grating. The guided modes of dielectric 
waveguides can be decomposed into plane waves. 3 This decomposition 
is particularly simple in the case of the modes of a slab waveguide. 
The following discussion is thus applied to this structure. Two plane 
waves are superimposed to form a traveling wave in the z direction-the 
direction of the waveguide axis-and a standing wave in the direction 
transverse to the z axis. The coupling coefficients for guided mode 
coupling and the radiation loss coefficients can be calculated by solving 
the plane wave scattering problem at the dielectric interface. 4 The 
geometry of the problem is shown in Fig. 1. For clarity of discussion, 
it was assumed that the incident plane wave approaches the interface 
at right angles. The actual plane waves making up the guided mode 
of the slab waveguide approach the interfaces at grazing angles. Figure 1 
is drawn with a sinusoi dally distorted core-cladding interface. In this 
case, the incident plane wave decomposes into a wave that continues 
to travel in the original direction after passing the interface and into 
a reflected wave plus a series of side lobes that are labeled by positive 
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Fig. 1 — A sinusoidally deformed dielectric interface functions as a phase grating. 
The figure shows the propagation vector of an incident plane wave (labeled 0) and 
reflected as well as transmitted plane waves of the higher-order grating lobes. 
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and negative integer numbers. These are the lobes of first, second, and 
higher order of the phase grating. If the incident plane wave make as 
grazing angle with the interface, the situation is essentially the same; 
the only difference being that some of the grating lobes form imaginary 
angles and are consequently evanescent instead of traveling waves. 
The zero-order lobe does not pass if the incident wave meets the inter- 
face at less than the critical angle for total internal reflection. Its 
angle is thus imaginary and only an evanescent wave exists on the far 
side of the interface while the incident wave is strongly reflected. The 
transmitted first-order grating lobes may both have imaginary angles. 
In this case, there are only scattered reflected waves in the core of the 
waveguide that can combine to a new guided mode provided that their 
angles correspond to one of the allowed directions for guided modes. 
If one of the transmitted first-order side lobes emerges on the far side 
of the interface with a real angle, it causes power to radiate into the 
space outside of the waveguide core. This radiation is lost to the guided 
wave and must be counted as power loss. 

For small amplitudes of the sinusoidal core-cladding interface irreg- 
ularity, the amplitudes of the grating lobes decrease rapidly with 
increasing grating order. The first-order grating lobe corresponds to 
the contribution of first-order perturbation theory and is indeed the 
only contribution of real interest if the core-cladding interface irreg- 
ularity has a small amplitude. With increasing amplitude of the interface 
irregularity, the higher-order grating lobes become increasingly im- 
portant. Their importance is enhanced by the fact that the angles of 
the second-order grating lobes may be real even when the transmitted 
first-order lobes both have only imaginary angles. This means that, to 
first-order of perturbation theory, no scattering loss exists. It is thus 
necessary to study the higher order grating responses in order to obtain 
information about scattering losses in case the first-order theory predicts 
no scattering loss at all. 

The study of these higher-order grating lobes and the derivation of 
power loss coefficients for the higher-order processes is the object of this 
paper. It is possible that the grating problem has been solved before 
to the accuracy that is attempted here. Because of the enormous volume 
of literature that exists on scattering problems, relevant papers may 
have escaped the author's attention. However, the application of the 
grating theory to the waveguide loss problem is probably new. 

We derive coupled equation systems for the amplitudes of the grating 
lobes and use these exact equations to obtain first- and second-order 
approximations in a straightforward way. Perturbation solutions of 
the exact equation system are particularly appropriate, since each 
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higher-order solution can be computed from the known lower-order 
solutions with no need to recompute the lower-order solutions each time 
the order of perturbation theory is increased by one. Formulae for the 
first- and second-order scattering loss process from a sinusoidal core 
cladding interface irregularity are derived. The results of this paper 
will be used in a companion paper 5 to calculate the loss penalty for 
intentional mode mixing in multimode waveguides. 

II. PLANE WAVE SCATTERING AT A SINUSOIDAL INTERFACE 

We consider the problem of a plane wave that impinges on the 
interface between two dielectric media. The interface is described by 
the function 

j{z) = a sin <£s. (2) 

If only first-order scattering is considered, more general shapes of /(z) can 
be synthesized by superposition of sinusoidal functions. For higher-order 
processes, mixing of the sinusoidal terms occurs so that the descrip- 
tion of scattering from more general interfaces becomes complicated. 
The incident plane wave is given by (the time dependence is e lut ) 

E v = Ae- iUix ^ iz) (3) 

H z = - &- Ae- iUi * +f>1,) (4) 

H M = ^-Ae~^ ix ^ iz \ (5) 

0}fl 

The coordinate system is shown in Fig. 1. It is assumed that no variation 
of either the field components or the material parameters exists in y 
direction so that we can symbolically write 

Ty " °- » 

The remaining three field components E x , E z , and H v vanish. The 
parameters k, and /3 ; are connected by the following equation 

« = «k a - tf) 1/2 . (7) 

The refractive index n Y belongs to the medium from which the plane 
wave approaches the interface, and k is the free space propagation 
constant. 

The reflected and scattered waves are expressed as superpositions of 
plane waves. We thus have in medium 1 
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E v = f BidY^-^ dp (8) 

J— to 

H, = -— f WW 1 '"'* dfi (9) 

H, = -— /"° o-.B(|8)e , ' t * lH '' ) <#, (10) 



with 



a = (n{k* - ff'\ (11) 
and similarly in medium 2 

E„ = f C03)e- , ' ( ' ,I+^I, dp (12) 

ff, = -— f* pCW* 1 "*** dp (13) 

h, = - C papy* 1 ""* m, (u) 



with 



- (# - fl') 1 ". (15) 



The boundary conditions at the dielectric interface require the 
tangential component of the electric field E„ and the tangential com- 
ponent of the magnetic field 

1 f 

H * = (! + jy H > + (i + f*p H * ( 16 ) 

to be continuous [/' is the z derivative of the interface function (2)]. 
The boundary conditions thus lead to the two equations: 

Am-'*'*'**' 1 + f* BW"-"'* dp = f C(0)e- i(fif+0l) dp (17) 

J -ao •'—go 

and 
(«, - fitf)AM- ttBt '+"* -["(« + 0/OBOy { * /_w dP 

J —ao 

= r (p- PnC(P)e- i(pf ' ei) dp. (18) 
In order to remove the z dependence from the equation, we multiply 
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with exp (—ifi'z) and integrate over z from — as to ». This procedure 
transforms the equations to the following form: 

AF(«, , ft - 00 

OO 

+ /"" B(0)F(-<r, - 00 d0 = f C&Wp, P ~ 0') #, (19) 
and 

+ f B<0)<?(-«r, 0, 00 <20 = f C(P)G( P , jS, /SO d0 (20) 

J -to "—no 

with 

F(,, |S - 00 - P B - ww+tf -^ > " dz, (21) 

J -co 

and 

G( v , 0, 00 = f" h - mz)Y iUni) + { ^' U ' dz. (22) 

J —00 

It is shown in the Appendix that F and G are related in the following 
way: 

G{ n , 0, |8') = *' + ^ = ^ F( v , 0-00- (23) 

Substitution of (2) into (21) yields 

Fin, - 0') = f" e --" — *v' ( ^'" c&. (24) 

The first exponential function under the integral sign can be expressed 
as a series in terms of Bessel functions with the help of the generating 
function of the Bessel functions. The remaining integration yields delta 
functions so that we obtain 

F(v, - 00 = 2ir £ JJM s tf - 0' + **]■ ( 25 ) 

» = — 00 

Substitution of (23) and (25) transforms the equation systems (19) 
and (20) to the form 

r-— oo 

= A £ J.M) 503. - 0'+W>), (26) 

|i=— 00 
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and 

.• = -00 I. &V P-y J 

= A £ ^^ ~ ^ /^a) 8(|3, - 0' + *). (27) 

v = — go *i 

The following abbreviations were used 

ft, = 0' + ^ (28) 

«r, = (# - #)* (29) 

and 

p, = (n 2 2 k 2 - ffl. (30) 

We know from the theory of phase gratings 2 that only discrete plane 
waves appear in the reflected and transmitted beams. This and the 
appearance of the delta functions on the right-hand side of (26) and (27) 
suggest that the solutions should be of the form 

OS 

B0) = E K 803 - ft - fit), (31) 

and 

C(0) = E C.803-& - tf,). (32) 

fi**— cO 

Substitution of (31) and (32) into (26) and (27) and comparison of the 
coefficients of the delta functions of equal arguments leads to two 
infinite equation systems for the unknown coefficients &„ and c„ . 

DO 

E {-b„ +r J f (a n+r a) + c n - r J.( Pn _ r a)\ = AJ^^a), (33) 

and 

V / ff " ~ "*(& + "*) J, w „> 

r-— oo k °"n + H 

Pn + y»(ft + «<ft) r , \ k\ - rufrpj ATf . f „ AS 

H c n -, r J,{p n - r a)? = — AJ n (*c,a). (34) 

Pn-v J K,- 

The equation system (33) and (34) is exact. An exact solution appears 
impossible to obtain. However, the equation system is very convenient 
for obtaining perturbation solutions of arbitrary order. It is also possible 
to obtain a solution that is exact in the limit <f> — > 0. 
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The reflected and transmitted fields follow directly from (8), (12), 
(31), and (32). In medium 1 we obtain for the reflected field 

B„ = £ b& il "*->>'\ (35) 

li — — to 

The field in medium 2 is 

E u = £ (vf""""**". (36) 

jt— — 06 

The parameters <r„ and p M are defined by (29) and (30). However, now 
we must use these equations with /3 M = £, : + m#- 

IH. SOLUTION FOR <£ — ► 

In the limit <f> — > 0, an exact solution of the equation systems (33) 
and (34) can be obtained. If = is assumed, we use the facts that 

<r v = k { (37) 

p„ = po (38) 

to write (33) and (34) in the form 

£ {-b n+ J r (Kia) + c n .J y ( PQ a)} = AJ n (KiO) (39) 

V = — CO 
BO 

X) {Kib n + y J v (KiO) + PoCn_,/,(^o)j = K t AJJM* ( 40 ) 

r = — 00 

It is now possible to eliminate c„ from the equations and obtain an 
equation system for b„ alone, 

f: b n+ jM = *-^^ ajm. (4D 

,rr« k» + Pa 



Similarly, we obtain by eliminating b t 



2& t 



£ c n -»J r (p a) = *' AJJi^a). (42) 

These equations can be solved with the help of the addition theorem for 
Bessel functions 

£ J n+ .(x)J,<3/)e M = UR)e in \ (43) 

p» — 09 

with 
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R = (x 2 + if - 2xy cos 8) u \ (44) 

With = 0, we obtain R = x — y so that we see immediately that 

&„ = Kj ~^- AJ&w) (45) 

Ki ~r po 

is the solution of (41). Since J-„(r) = J v (— x) (for integer values of v) 
it is also apparent that 

2k 

c„ = — -r~ AJ^ - po)a] (46) 

Ki T Po 

is the solution of (42). 

The solutions (45) and (46) are exact for <j> = 0. One might expect 
that $ = describes a plane dielectric interface so that no side lobes 
should be expected. Even though it is true that all the sidelobes coincide 
for <j> = 0, the solutions (45) and (46) do hold approximately even if 
t}> ^ 0. The sinusoidal shape of the interface is apparently built into 
the equation system (41) and (42) even though $ does not appear 
explicitly. The solutions (45) and (46) are approximations that hold 
if hj> <K 1. These solutions show that the amplitudes of the side lobes 
are proportional to Bessel functions. This result is well known from the 
theory of phase gratings. 2 

IV. PERTURBATION SOLUTIONS 

For our purposes, the solution for hf> « 1 is of little use. Therefore, 
we proceed to derive approximate solutions that hold for 

ha « 1. (47) 

We use the following approximations for the Bessel functions of small 
argument 

J Q (x) = 1 - f (48) 

Ux) = -J-&) = |(l-|-) (49) 

J 2 (x) = J_ 2 (x) = |- (50) 

In addition, we assume that b„ and c„ are zero-order terms, b±i and c ±l 
are of first order, and b ±2 and c i2 are of second order. 

By neglecting all but zero-order terms in (33) and (34), we obtain 
to zero order of approximation for n = 0, 
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— b + Co = A 
ffobo + PcA> = K iA 
Taking n = 1, we obtain to first order 



— &! + c, =^ (k, A — 0" &o — poCo) 



(51) 



- [pi + ^03, + *)]«>} 



(52) 



The corresponding equation system for b-i and c. x is obtained by 
using n = — 1. It has exactly the same form as the equation system (52) 
and is obtained by replacing bj with — 6_i , c x with — c_i and with —<p. 
Finally, we obtain the second-order approximation by setting n = 2 
and keeping only terms up to second order 

-6a + C 2 = ^ (K t A + o-o&o - PoCo) - ^ (<Ti&i + PiCi) 

<r 2 b 2 + p 2 c 2 - I Mr? - 2*0,-)ii - ffoW + 2*08, + 2<p)]6 
- M + 2003,- + 2*)]co) 

+ § {[*z + *(ft + 2*)]b, - [p* + *(& + 2*)]c) 

The equations for b- 2 and c_ 2 are obtained by replacing b 2 with b_ 2 , 
c 2 with c- 2 , &i with — &_, ,c, with — c_, , and finally with -0 in (53). 

It is immediately apparent that each order of approximation follows 
from the preceding order. We can thus solve all the equations (51), 
(52), and (53) in succession. Each time we need solve only two equations 
with two unknowns. The result of the previous approximation is then 
used to obtain the next higher order of approximation from the next 
equation system. 

The solutions of these equations are listed below. 



Ki + Pa 

m *<(«.- ~ go) aA 
ffi + Pi 



Co = 



2*« 



*, + Pa 
C x = &i 



(54) 
(55) 
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h - 2-A„ Ki ~ Po L _L „ .9/ ^ «,- + Po T 

4 <r 2 -f- p 2 L o"i + PiJ 

I Po " IT, + 2(«, - p ) «L±Jo] 

L ffi + piJ 



C 2 = —r~ k £ 



*. — Po 

0"2 + p2 



(56) 



The coefficients &_! and c_j are obtained from (55) by changing the sign 
in front of the terms and changing the subscripts 1 to — 1 on the right- 
hand side of the equation. The signs of 6_ 2 and c_ 2 are the same as those 
of the coefficients in (56). We obtain these coefficients by changing 
the signs of the subscripts on the right-hand side of the equations. For 
ha « 1 and k<f> « 1 the equations (45) and (46) can be shown to be 
identical with (54), (55), and (56). 

The amplitudes may belong to plane traveling waves or to evanescent 
waves. Whether a wave is of the propagating or evanescent type depends 
on whether the parameters <r„ and p A are real or imaginary. For real 
values we obtain traveling waves while imaginary values indicate that 
the field decays exponentially with increasing distance from the interface 
indicating an evanescent wave. The propagation constant in z direction, 
j3 M , is obtained from (28) by replacing j8' with & . We thus have for 
traveling as well as for evanescent waves 

ft, = ^ + W> M = 0, ±1, ±2, • ■ ■ . (57) 

V. CALCULATION OP SLAB WAVEGUIDE LOSSES 

Since the guided modes of the slab waveguide can be expressed as the 
superposition of two plane waves, whose propagation vectors form 
equal but opposite angles with the z axis, we can use our present results 
immediately to calculate the radiation losses suffered by the guided 
slab waveguide modes. 4 Our calculation applies to TE modes. However, 
for slight index differences the losses of TE modes and TM modes are 
nearly identical. The slab waveguide geometry is shown in Fig. 2, 

The radiation losses of slab waveguide modes have a somewhat 
complicated dependence on either frequency or slab width, since the 
interference of the waves scattered at one of the two dielectric interfaces 
with the radiation from the other interface-and also the interference 
with radiation that is reflected at the opposite interface-has to be taken 
into account. However, these interference effects cause only fluctuations 
about an average value. If we content ourselves with establishing only 
the average loss value, disregarding the fluctuations, the description of 
radiation losses is greatly simplified. We also gain the advantage of 
obtaining simpler mathematical expressions. In the spirit of this sim- 
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Fig. 2— Cross section of the slab waveguide. 

plification, we consider all scattered power as lost- with the exception of 
those waves that are captured in the core-disregarding incomplete 
reflection from the other dielectric interface and interference with the 
directly scattered radiation. The amount of power scattered per unit 
length of the waveguide is 

S, = | | CBi X HTJ. I + h I (E* X WSk V (58) 

The electric and magnetic fields in this expression are only the scattered, 
untrapped part of the field, exclusive of the incident field. The sub- 
scripts 1 and 2 refer to the fields in mediums 1 and 2. With the assump- 
tion that only one side lobe is instrumental in dissipating power by 
radiation we can write 

S x = ^-(^ | &., | 2 + p-, | ol, ft. (59) 

zw/x 

The subscript v assumes the values 1 and 2 for first- and second-order 
light scattering. The negative values must be used (v is now assumed to 
be positive) because o-_, and p-„ must be real since evanescent waves 
do not carry power. 

The power attenuation is 

- = sb- (60) 

S, is the power pe runit length (unit area in the three dimensional case) 
that is carried by the plane wave in z direction. The total power carried 
by this one plane wave component in z direction inside of the waveguide 
core of width 2d is thus 2S,d. The ratio of the power lost (per unit 
length along the waveguide axis) divided by the power carried by the 
wave is the power loss per unit length. Actually, two plane waves are 
needed to describe the guided mode in the slab waveguide. However, 
thus far we have considered only the scattering loss from one of the 
two interfaces. It is sufficient to consider that each of the two plane 
wave components scatters from one interface. Inclusion of the scattering 
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from both interfaces introduces a factor of two in the numerator of (60). 
However, another factor of two is introduced in the denominator by 
adding the power of the other plane wave to the total power carried by 
the guided mode. The expression (60) thus holds for the scattering 
loss of the guided modes provided that both interfaces contribute an 
equal amount to the power loss. Using the following expression for the 
power flow density in ^-direction, 

S - ~ 4 I A I"' (61) 

we obtain the general expression for scattering losses from a dielectric 
slab waveguide with sinusoidally deformed core-cladding interfaces 



1 

2&< 



~ ? . | A fig (<r- | &-, I" + P-, | C_„ | 2 ). (62) 



The index v indicates the order of the scattering process. Taking v = 1, 
we obtain with the help of (55) the scattering loss contribution from 
the first-order grating lobes 

2/ 2 2\i 2 

_ K,(n 1 - n 2 )k a 
011 ~ orT — T vj a • ( 63 ) 

We used the fact that k< is real, while p is imaginary, so that we have 
I *i + p T = k\ — pi = («i — n\)k 2 . 

The loss contribution of the second-order side lobes follows from (56) 
and (62) with v = 2. 

" - mitr+%i ^ + '-*- + *' - 2 "- n (64) 

The 7 parameters are defined by the equations 

7< = i Po = (fi - n 2 2 k 2 ) l/2 (65) 

and 

T_, = ip. t = [flfc - 0) 2 - nlk 2 ] l/ \ (66) 

We have tacitly assumed that the first-order side lobe belongs to an 
evanescent wave if the loss contribution of the second-order side lobe 
is being considered. The parameters y { and y- v of (65) and (66) are 
consequently assumed to be real quantities. If the first-order side lobe 
propagates in medium 2 as a traveling wave, the second-order side lobe 
also gives rise to a traveling wave. However, since the contribution 
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from the grating lobe of first order is much stronger than that from the 
second-order lobe, we can neglect the loss contribution (64). If first- 
order scattering does occur, it is the predominant effect. Only if first- 
order scattering contributes only an evanescent wave, and does not 
cause radiation loss, must the second-order loss process (64) be con- 
sidered. 

VI. DISCUSSION 

We have calculated the loss contributions that result from the grating 
lobes of first and second order when a plane wave impinges on the 
sinusoidally deformed interface between two different dielectric media. 
Not all grating orders belong to traveling waves. Some grating orders 
cause evanescent waves in medium 2 and guided waves in the core so 
that they do not contribute to radiation loss of a guided wave in medium 
1. The first-order radiation loss coefficient (63) is proportional to (ah) 2 
while the second-order radiation loss coefficient (64) is proportional 
to (afc) 4 . If both processes are effective simultaneously, the lower order 
process is dominant. Whether the second-order process is the only cause 
of radiation loss, or whether both first- and second-order processes are 
acting simultaneously depends on the magnitude of the mechanical 
frequency of the sinusoidal interface distortion. If 7-1 of (66) is real, 
the first-order side lobe belongs to an evanescent wave, and only the 
second-order side lobe is causing radiation losses. When y^ is imaginary, 
both first- and second-order side lobes carry away real power. Equation 
(64) is not applicable in this case since it was derived under the assump- 
tion that (65) and (66) are both real. However, if first-order radiation 
losses are possible, the second-order loss coefficient gives only a small 
contribution to the total radiation loss. 

It is interesting to compare the result of our present theory with 
earlier results obtained from a modal analysis of the slab waveguide 
problem. From equation (79) of Ref. 1, we obtain in our present notation 

_ ak*(7\\ — rily* |~ p-i cos 3 a^d 

/ 1 \ Lp-i cos 2 <!-$ + cr-i sin 2 tr-jd 

J,d(j ' 



4j8,-d 1 + 



7 



id) 



■ p., sin 2 <r. x d 1 ffi? v 

plj sm ff-id + tr_i COS cr^dj 



Equation (67) is modified for the case that both interfaces are sinusoi- 
dally distorted, but with a random phase relationship between the two 
sinusoidal functions. The radiation loss (67) is considered to be an 
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ensemble average over slab waveguides with all possible phase relation- 
ships between the sinusoidal distortions of the two interfaces. In 
addition, we used the expression 

z 

cos * «' d = („; - V (68) 

which represents the eigenvalue equation of the even TE modes of 
the slab waveguide. 

It was shown in equations (39) and (40) of Ref. 6 that the average 
value of the expression in brackets of eq. (67) is given by 

2/(<r_ r + p_,). (69) 

Combining (69) with (67) makes it apparent that the average value of 
the slab waveguide radiation loss is identical to the loss coefficient (63) 
that was derived from the plane wave model. The only remaining 
difference can be explained as stemming from the fact that the actual 
width 2d of the slab must be replaced by the effective slab width 

which is caused by the exponential field tail reaching out into the 
cladding. 

The plane wave model used in this paper does not include the inter- 
ference effects of waves originating from the two interfaces, and from 
the reflection of the scattered light from the opposite interface. But the 
averaging processes that were involved in converting the precise 
scattering loss coefficient of the slab waveguide theory into the loss 
coefficient derived from the plane wave model may be expected to be 
effective in a real waveguide. If the phase of the sinusoidal interface 
distortion varies slowly and randomly along the waveguide, the phase 
average that is already incorporated in (67) may actually occur. The 
average over the expression in brackets in (67) would occur either with 
randomly changing slab half width d, or with randomly varying mechan- 
ical frequency <f> of the sinusoidal interface changes. The loss formula (63) 
is much simpler than (67), and is of actual practical value for the 
indicated reasons. The same averaging process involved in the first-order 
loss coefficient (63) is also implicit in the second-order loss coefficient 
(64). 

Second-order scattering couples two modes v and n, if their propagation 
constants satisfy the condition 

I A - A | - 20. (71) 
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Mode v is a guided mode, while mode p can be either a guided or a 
radiation mode. Equation (71) replaces eq. (1) for first-order scattering. 
Higher-order processes have similar coupling laws with the order of 
the process multiplying the mechanical frequency <f>. 

Second-order scattering losses are of importance for intentionally 
coupled multimode operation. Mode coupling reduces the pulse delay 
distortion that is caused by the different group velocities of the modes. 
It is possible to design the core-cladding interface irregularities in such a 
way that all guided modes (with the exception of the last) are coupled 
to each other without coupling to the continuous spectrum of radiation 
modes by first-order processes. 7 However, radiation losses via second- 
and higher-order processes are still possible. The discussion of second- 
order losses for intentionally coupled multimode operation is the subject 
of a companion paper. 5 

APPENDIX 

Proof of the Relation (23) 

We begin by using the fact that a sinusoidal function with an infinite 
argument can be regarded as zero. This assertion is the basis for the 
following definition of the delta function: 

5(,0 = lim± S -^^ (72) 

A^x, IT X 

vanishes everywhere except at the point x = 0. The singularity of the 
delta function at x = is caused by the appearance of x in the denomi- 
nator. Without this denominator, we are justified to define 

lim sin xA = 0. (73) 

A-to 

Using eq. (73), we write the following identity 

- lim - 2isin [ v f(A) + (fi - 0*)A] 

= lim T Ae-*'"*'^-'' 1 " dz 

A~m J-A O.Z 

- i[vf + - ffi-wv-™ dZm (74) 



-L 



It was assumed that f(-z) = —f(z). From the last line of (74), we 
obtain immediately 
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or 

f h-/3/^)]e-'' [ " /(l, + ( ^' )l1 ^ 

■r — oo 

= £ + w - f) r ,_,„„.„„_,.,., ffe 

7 J — no 

From the definitions (21) and (22) it follows that (75) is identical 
with (23). 
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